Skyrmion Strings and Anomalous Hall Effect in Double Exchange Systems. 
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We perform Monte Carlo simulations to obtain quantitative results for the anomalous Hall re- 
sistance, Ra, observed in colossal magnetoresistance manganites. Ra arises from the interaction 
between the spin magnetization and topological defects via spin-orbit coupling. We study these 
defects and how they are affected by the spin-orbit coupling within the framework of the double 
exchange model. The obtained anomalous Hall resistance is, in sign, order of magnitude and shape, 
in agreement with experimental data. 
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Doped perovskite manganites have attracted much 
attention lately, since they show the so called colos- 
sal magnetoresistance Q. These materials undergo a 
ferromagnetic-paramagnetic transition accompanied by a 
metal-insulator transition. The Double-Exchange (DE) 
mechanism Q plays a major role to explain this mag- 
netic transition. In the DE picture, the carriers moving 
through the lattice are strongly ferromagnetically cou- 
pled to the Mn core spins and this produces a modula- 
tion of the hopping amplitude between neighboring Mn 
ions. 

The Hall resistivity pu in ferromagnets has two con- 
tributions [^, one proportional to the magnetic field H, 
and the other to the spin magnetization M: pH=RoH + 
RaM. Ro and Ra denote the ordinary (O) and the 
anomalous (A) Hall resistances (HR). The existence of 
Ra requires a coupling of orbital motion of electrons to 
M, and the AHE is usually explained in terms of skew 
scattering due to spin-orbit (s-o) interaction [Q. 

Several groups |^-|| have measured pn of the doped 
Mn oxide LaQ -jCa^^^MnO^ for different temperatures 
(T). These experiments found that Rq and Ra have 
opposite sign, that Ra is much bigger than i?o, that Ra 
peaks at a T above Tc, and decreases slowly at higher T. 
These effects can not be explained with the conventional 
skew scattering theory Recently, it was proposed 
[^,0that in Mn oxides Ra arises from the interaction 
of M with non trivial spin textures (topological charges) 
via s-o coupling. The number of topological charges in 
the three-dimensional (3D) Heisenberg model increases 
exponentially with T and Ye et al. extrapolated 
these results to DE materials obtaining a rapid increase 
of Ra at low T. For T > although in 3D there is 
not a theory of topological defects. Ye et al. were able to 
estimate the overall shape of Ra- 

In this work we perform Monte Carlo simulations in 
order to obtain a quantitative form of Ra in DE sys- 
tems. We compute the number of topological defects as 
a function of T and, by introducing s-o interaction, we 
couple the defects orientation with M. In this way we 
obtain an AHR which has the same sign, shape and order 
of magnitude as found experimentally 

Hamiltonian. The electronic and magnetic properties 



of the Mn oxides are described by the DE Hamiltonian, 

-^Jh^ d^^^(T^^„,di^^'Si + Hso + , (1) 

here df^ creates an electron at site i and with spin a, t 
is the hopping amplitude between nearest-neighbor sites, 
Jh is the Hund's rule coupling energy. Si is the core 
spin at site i and Hso is the spin-orbit (s-o) interaction. 
In the tight binding approximation, the effect of H is 
to modify the hopping matrix element by introducing 
the phase a^A^g, where a is the lattice parameter, j — 
i + 5 and A is the vector potential corresponding to H. 
The last term in the Hamiltonian is the Zeeman coupling 
Hz = 5MbH^j Si. We assume that the Mn ions form a 
perfect cubic lattice. 

In the limit of infinite J^f, the electron spin at site i 
should be parallel to Si, and H becomes 



H = 



't y~cos 



^g<0(«j)/2+ajfeyl,j) + 



Hso+Hz. (2) 



Now cf creates an electron at site i with spin parallel 
to Si, mi=Si/iS' and cos9ij—mi-mj. The first term in 
Eq.(2) describes the motion of electrons in a background 
of core spins. The electron hopping is affected by the 
nearest neighbor spin overlap, cos(6'ij/2), being the ki- 
netic energy minimum when all core spins are parallel. 
This is the DE mechanism for the existence of a ferromag- 
netic metallic ground state in Mn doped oxides. When 
the orientation of the electron spin is moved around a 
close loop the quantum system picks up a Berry's phase 
|]l3| proportional to the solid angle enclosed by the tip of 
the spin on the unit sphere. 4>{i,j) is the Berry's phase 
defined mathematically as the solid angle subtended by 
the unit vectors nii, irij and z on the unit sphere. 

The phase (j>{i,j) in the Hamiltonian affects the motion 
of electrons in the same way as does the phase arising 
from a physical magnetic field. 4>{i,j) is related with in- 
ternal gauge fields generated by non trivial spin textures, 
which appears in the system when increasing T. In ab- 
sence of s-o coupling the phases 4){i,j) are random and 
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the net internal gauge field is zero. In the presence of s-o 
coupling there is a privilege orientation for the spin tex- 
tures and a non zero average internal gauge field appear. 

In the following we study, as a function of T, the ap- 
pearance of topological defects in the DE model. Once 
the spin textures are characterized, we study the effect 
that the s-o coupling has on the defects and we analyze 
their contribution to the Hall effect. 

Topological defects in the DE model. The temperature 
induces fluctuations in the ferromagnetic state of the core 
spins and at T > Tc the system becomes paramagnetic. 
Typically Tc ~ 300i^ and this T is much smaller that the 
electron Fermi energy, t ~ O.leT^, therefore we consider 
that the conduction electrons temperature is zero. 

In first order in the electron wave function the temper- 
ature dependence of the magnetic properties is described 
by the classical action |14|, 



S^l3t{c-^ 



^ j)/2 



(3) 



where ( )o means expectation value at T = 0. 

The system described by Eq.(3) is expected to behave 
with T similarly to the Heisenberg model. In particu- 
lar, at low r, we expect the occurrence of point singu- 
larities as those occurring in the 3D Heisenberg model 
These points singularities are called topologically 
stable because no local fluctuations in a uniform system 
can produce them. The defects are classiflcd by the 
topological charge, Q, which represents the number of 
times and the sense in which spins on a closed surface 
surrounding the defect cover the surface of a unit sphere 
in spin space. An example of topological defects with 
Q = 1) occurs at a position from where all the 

spins point radially outward (inward). 

In order to locate the topological charges in the lattice 
model we follow the prescription of Berg and Liisher . 
For each unit cube of the lattice we divide each of its six 
faces into two triangles. The three unit-normalized spins 
at the corners of a triangle I define a signed area Ai on 
the unit sphere. The topological charge enclosed by the 
unit cube is given by 



1 

q^—Ta 



(4) 



1=1 



This definition of Q ensures that, when using periodic 
boundary conditions, the total topological charge in the 
system is zero. Q in each cube is an integer number, and 
the magnitude of nonzero charges is almost always equal 
to unity. Only at rather high temperatures a few defects 
with IQl > 2 are found. 

Q can be interpreted as the number of magnetic 
monopoles enclosed by the unit cube. The quantities 
Ai(f>o/4:TT represent the magnetic flux piercing the trian- 
gle I. Here (po = hc/e is the magnetic flux quantum. 



Following this analogy we assign at each point of the lat- 
tice, i, a three dimensional internal magnetic fleld b^. 
The phase 4'ihj) can be written in the form, 



e 

10 ' 



he 



(5) 



where b=V x a. Here it is clear that the b associated 
with a spin texture affects the motion of an electron just 
as does an external magnetic field. 

In a system with a uniform magnetization at the sur- 
face, the interaction between a positive and a negative 
defect is finite and in the continuous Heisenberg model 
increases linearly with separation. Jl6t Defects with op- 
posite charges are closely bound in pairs at low T. These 
pairs of defects are Skyrmion strings (dipoles) 

which begin at a monopole (Q^l) and end at an an- 
timonopole {Q=-l). The Skyrmions are characterized by 
a dipole P joining Q—-1, and Q=+l, 



(6) 



By performing Monte Carlo simulations on the classi- 
cal variables m; , we have studied the dependence on 
T of the number of defects in the system. With the defi- 
nition of Q, the number of positive and negative defects is 
the same and the important quantity is the average defect 
pair density (n) . In Fig. ^ we plot (n) as a function of T 
for a DE system of size 16x16x16. We have checked that 
our results arc free of finite size effects. In the same figure 
we plot m(T). In the DE system Tc ~ 1.06 t(c+Cj)o. The 
results correspond to gfibH=0 and gfibH=0.067Tc {H ^ 
10 Tesla). 
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FIG. 1. Plot of the logarithm of the average defect pair 
density as a function of the inverse temperature for a system 
of size 16x16x16. The dashed lines are linear fits to the data 
below Tc- Their slope is the core energy Ec — 7.05Tc. 

We have studied correlation between defects. For 
T < Tc the Q=+l and Q=-l defects are strongly cou- 
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pled forming Skyrmions. The Skyrmions are very dilute, 
almost independent, and their density can be fitted to 



in) 



(7) 



Here Ef. and a are, respectively, the core energy and the 
entropy of the Skyrmions. Numerically we have obtained 
Ec = T.OSTc- This value is slightly smaller than the value 
obtained for the Heisenberg model Q E^"'' = 8.7 T^. 
At T > Tc the number of defects increases sharply and 
it becomes very difficult to pair up defects with opposite 
charges in an unambiguous way. The core energy practi- 
cally does not depend on H, this result is in agreement 
with the obtained for pure two-dimensional Skyrmions 
^ . The entropy a is related with the degeneracy of the 
Skyrmions in the orientation of P and with thermally 
activated twist and dilatation soft modes. We have ob- 
tained that for H=0, a=5l. Note that a""'' - 320 @. 
For finite H the six (1,0,0) directions keep being degen- 
erated, but some of the twist soft modes becomes more 
gaped and the value of the entropy is reduced. 

By performing calculations in a system constrained to 
get a unique Skyrmion with a given dipole P, we have 
calculated the energy of isolated Skyrmions. In this way 
we have checked that in the absence of s-o coupling the 
energy of the Skyrmion only depend on the absolute 
value of P and not on its orientation. We have obtained 
that the core energy of a Skyrmion with P=(1,0,0) is 



E 



(1,0,0) 



=7.05Tc, with P=(1,1,0) is E, 



(1,1,0) 



=8.65 Tc and 



for P=(l,l,l) is e!^''^'^'' =:10.11 T^. The dependence of 
the energy on P confirm the strong confinement energy of 
the topological defect . The increase of the Skyrmion 
energy with P, implies that for T < Tc the only relevant 
pairs are those with P—1. 

In the case of zero s-o coupling, the Skyrmion core 
energy does not depend on the direction of P and the 
average internal gauge magnetic field is zero, (b)=0. The 
s-o coupling privileges an orientation of the Skyrmions 
and it results in a finite value of (b). 

Spin-orbit interaction. The s-o coupling has two con- 
tributions, the interaction of the spin of the carriers with 
the ion electric field, and the magnetic interaction be- 
tween the carriers with the core spins Si . In the Jh — > oo 
limit, both contributions have the same form and. 



Hso — ^soT^^S nii ■ hi 



(8) 



where x is the carriers concentration. For Aso ^ and 
m ^ 0, the Skyrmions are preferentially oriented. If 
H II z, m II z and a net z component of (b) results. 

We have included the s-o interaction in the Monte 
Carlo simulation and we have obtained that (62) is linear 
with \so in all range of T . In Fig. ^ we plot for a sys- 
tem of size 12 X 12 X 12, 4'z = {bz)a'^ / 4>o as a function of T. 
For T < Tr our results are free of finite size effects. For 



T > Tc the increase of the system size does not change 
the overall shape. 
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FIG. 2. Plot of (j>z as a function of T for a system of size 
12 X 12 X 12 and gfibH = 0.067Tc. The interaction of the spin 
magnetization and topological defects via spin-orbit coupling 
has been included. The dashed line corresponds to an inde- 
pendent Skyrmions picture (see Eq.|ll|). 

It is interesting to analyze the results for (j)^ as a 
function of the density of defects (n). The main effect 
of the s-o coupling is to privilege the appearance of 
Skyrmions polarized parallel to m, (n^), with respect to 
Skyrmions with P antiparallel to m, This asym- 

metry produce a z-component of the internal gauge field, 
(&z}=0o/fl^(("-) ^ ("■+))■ Assuming that the Skyrmions 
are independent. 



(n_)-(ri+) = -i(n)sinh(eo/r) , 



where 



£0 



^soX ^ , . 
— n — "'■(-' ) 



(9) 



(10) 



is the s-o energy interaction of a Skyrmion with P = 
(0, 0, 1) and rh{T) is the spin polarization inside the 
Skyrmion. Finally, within the independent Skyrmions 
picture, after linearizing Eq.|^, 



lind 



jn) Xso 
6 T 



X rh{T) 



(11) 



'm{T) is expected to be a function of m{T). By comparing 
the expression for with the Monte Carlo results. Fig. 
|, we obtain that for T < T^ m{T)^m{T)/5. That 
means that the spin polarization inside the Skyrmions 
is five times smaller that the average spin polarization. 
It is interesting to note how the independent Skyrmions 
picture describes not just the low temperature regimen, 
but also the trends of (b^) at T > Tc. 

The sign of (f)^ depends on the sign of Aso- Physically 
we expect that the motion of the electrons leads to an 
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internal gauge field which acts to cancel the applied field, 
and therefore a negative sign for (6^) is expected. 

Anomalous Hall effect. The Hall resistivity can be 
written as pn = {H + {b))/nec. Comparing this expres- 
sion with the definition of Ra we obtain 



Ro gfJ-b m 



(12) 



By using the appropriate parameters, we plot in Fig. ^ 
Ra/Ro for Lao.TCao.sMnO^. We take - 2701^ [|||; 
we use Xso=5K @ and gfitH = OMlTc. 

Fig. 3 is the main result of our calculation. We obtain 
an AHR which i) is negative, ii) increases exponentially 
and becomes evident at T ~ Tc/2, iii)at temperatures 
close to Tc, is around 20 times bigger than the ordinary 
Hall resistance, and iv) has a maximum at temperatures 
slightly higher (~ 30K) than T^. The subsequent de- 
crease is due to thermal fluctuations that destroy the 
directional order of the Skyrmions. If we take into ac- 
count the fact that the conductance at T > Tc is mainly 
due to polaron hopping, an extra factor of is ex- 

pected [1^ leading to a steeper decrease. Our results are 
in good agreement with the data obtained experimentally 
We have checked that these results do not depend 
significantly on the magnetic field applied. 

In addition, we have found, by diagonalizing the elec- 
tron Hamiltonian, that there is not significant electronic 
charge associated to the Skyrmions |Q . This is in con- 
trast to the quantum Hall ferromagnetic systems where 
the topological and electrical charge are equivalent [ p^ . 
The Skyrmions we have studied appear when increas- 
ing T but we have found that they can also appear at 
T = when spin defects (in particular, antiferromagnetic 
islands) are present in the system [pO[. 
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FIG. 3. Plot of Ra/Ro as a function of temperature. This 
is our main result. We recover the sign, order of magnitude 
and shape of the measured anomalous Hall resistance. 



sistance in colossal magnetoresistance manganites. We 
have obtained, as a function of temperature, the aver- 
age value of Skyrmions for the double-exchange model. 
By introducing a spin-orbit interaction, we obtain a cou- 
pling between the Skyrmions orientation and the magne- 
tization which results in the appearance of an anomalous 
Hall resistance. The sign, order of magnitude and shape 
of the obtained anomalous Hall resistance are in agree- 
ment with the experimental information. 
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